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ABSTRACT 



Lightcurves have long been used to study stellar activity and have more recently become a major tool in the field of 
exoplanet research. We discuss the various ways in which stellar activity can influence transit lightcurves, and study 
the effects using the outstanding photometric data of the CoRoT-2 exoplanet system. We report a relation between 
the 'global' lightcurve and the transit profiles which turn out to be shallower during high spot coverage on the stellar 
surface. Furthermore, our analysis reveals a color dependence of the transit lightcurve compatible with a wavelength- 
dependent limb darkening law as observed on the Sun. Taking into account activity-related effects, we re-determine the 
orbit inclination and planetary radius and find the planet to be ~ 3% larger than reported previously. Our findings 
also show that exoplanet research cannot generally ignore the effects of stellar activity. 
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1. Introduction 

The brightness distribution on the surface of active stars 
is both spatially inhomogcncous and temporally variable. 
The state and evolution of the stellar surface structures 
can be traced by the rotational and secular modulation of 
the observed photometric lightcurve. In the field of planet 
research, lightcurves including planetary transits are of par- 
ticular interest, since they hold a wealth of information on 
both the planet and its host star. 

The outstanding quality of the space- based photome- 
try p rovided by the CoRoT mission (e.g. lAuvergne et al.l 
I2OO90 provides stellar lightcurves with an unprecedented 
precision, temporal cadence and coverage. While primar- 
ily designed as a planet finder, the CoRoT data are 
also extremely interesting i n the context of stellar ac- 
tivity. Recently, iLanza et al.l (|2009l ) demonstrated the in- 
formation content to be extracted from such lightcurves 
in the specific case of CoRoT-2a. This star is solar-like 
in mass and radius, but rotates fa ster at a s peed of 
vsm{i) = 11.85 ± 0.50 km/s (Bouch v et al.ll2008[ ). Its ro- 
tation period of « 4.52 days was deduced from slowly evolv- 
ing active regions, which dominate the photometric varia- 
tions. Thus, CoRoT-2a is a very active star by all stan- 
dards. Even more remarkably, C oRoT-2a is orbited by a 
giant planet (|Alonso et al.l [200^ . which basically acts as 
a shutter scanning the surface of CoRoT-2a along a well 
defined latitudinal band. 

The transiting planetary companion provides a key 
to understanding the surface structure of its host star. 
While previous an alyses have either ignore d the transits 
(ILanza et al.ll2009 ] or the 'global' lightcurve (j Wolter et all 
I2009D . we show that there is a relation between the transit 
shape and the global lightcurve which cannot generally be 
neglected in extrasolar planet research. 



2. Observations and data reduction 



lAlonso et al.l ((lOOl) discovered the planet CoRoT-2b using 
the photometric CoRoT data (see Table [T]). Its host star 
has a spectral type of G7V with an optical (stellar) com- 
panion too c l ose to be resolved by CoRoT. According to 
lAlonso et al.l (|2008t) . this secondary contributes a constant 
(5.6 ±0.3)% of the total CoRoT-nieasured fiux. CoRoT-2b's 
orbital period of « 1.74 days is about one third of CoRoT- 
2a's rotation period, and the almost continuous CoRoT 
data span 142 days, sampling about 30 stellar rotations 
and more than 80 transits. The lightcurve shows clear ev- 
idence of strong activity: There is substantial modulation 
of the shape on time scales of several days, and the tran- 
sit profiles are considera bly deformed as a consequence of 
surface inhomogeneities (jWolter et al.ll2009l ). 

Our data reduction starts with the results provided by 
the CoRoT N2 pipehne (N2_VER 1.2). CoRoT provides three 
band photometry (nominally red, green, and blue), which 
we extend by a virtual fourth band resulting from the com- 
bination (addition) of the other bands. This 'white' band 
is, henceforth, treated as an independent channel, and our 
analysis will mainly refer to this band. It provides the high- 
est count rates and, more importantly, is less susceptible to 
instrumental effects such as long term trends and 'jumps' 
present in the individual color channels. 

In all bands we reject those data points flagged as 'bad' 
by the standard CoRoT pipeline (mostly related to the 
South Atlantic anomaly) . The last step leaves obvious out- 
liers in the light curves. To remove them, we estimate the 
standard deviation of the data point distribution in short 
(« 3000 s) slices and reject the points more than 3a off a 
(local) linear model. Inevitably, we also remove a fraction 
of physical data (statistical outliers) in this step, but we 
estimate that loss to be less than a percent of the total 
number of data points, which we consider acceptable. 

In all bands the white we find photometric discontinu- 
ities (jumps) which are caused by particle impact on the 
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Table 1. Stellar/planetary parameters of CoRoT-2a/b. 



Star ^ 


Value zt Error 


Ref._^ 


P. 


(4.522 ± 0.024) d 


L09 


SpBctrEil typ6 


G7V 


BOS 




Value zb Error 


xxer. 




(1.7429964 ± 0.0000017) d 


AOS 


n [BJD] 


(2454237.53362 ± 0.00014) d 


AOS 


i 


(87.84 ± 0.10)° 


AOS 


Rp/ Rs 


(0.1667 ± 0.0006) 


AOS 


a/Rs 


(6.70 ± 0.03) 


AOS 


Ua, Ub 


(0.41 ±0.03), (0.06 ±0.03) 


AOS 



" Ps - stellar rotation period. 

' taken from iLanza et al.1 (|2009l ') [L09], lAlonso et all l|2008f) 
[AOS], or lBouchv et al.l l|2008l ) [BOS] 

Pp - orbital period, Tc - central time of first transit, i - or- 
bital inclination, Rp,Ra - planetary and stellar radii, a - semi 
major axis of planetary orbit, Ua,Ub - linear and quadratic limb 
darkening coefficients. 



CoRoT detector. In the case of CoRoT-2a the jumps are of 
minor ampUtude compared to the overall count rate level, 
and we correct them by adjusting the part of the lightcurve 
following the jump to the preceding level. 

Finally, we correct the CoRoT photometry for system- 
atic, instrumental trends visible in all bands but white. In 
order to approximate the instrumental trend, we fit the 
(entire) lightcurve with a second order polynomial, q, and 
apply the equation 

^corrA — ^o.i ' • (1) 

Here, Co,i is the i-th observed data point and qi the associ- 
ated value of the best fit second order polynomial, c repre- 
sents the mean of all observed count rates in the band, and 
Ccorr,i the Corrected photometry. 

The resulting lightcurve still shows a periodic signal 
clearly related to the orbital motion of the CoRoT satel- 
lite. Again, this is a minor effect in the white band, and we 
neglect this in the context of the following analysis. 

In a last step we subtract 5.6% of the median lightcurve 
level to account for the companion contribution. Note that 
we use the same rule for all bands which can only serve 
as an approximation as lAlonso et al.l ()2008f ) point out that 
the companion has a later type (probably K or M) and, 
therefore, a different spectrum than CoRoT-2a. 

3. Analysis 

3.1. Transit profiles and stellar activity 

A planet crossing the stellar disk imprints a chara cteristic 
tran sit feature on t he lig htcurve of the star (e.g. Po nt et aTl 
[20071 IWolter et aLll2009t ). The exact profile is determined 
by planetary parameters as well as the structure of the stel- 
lar surface. A model which describes the transit profile must 
account for both. One of the key parameters of the sur- 
face model is the limb darkening law. The presence of limb 
darkening seriously complicates transit modeling, because 
it can considerably affect the transit profile, while it is hard 



to recover i ts characteristics from lightcurve analyses (e.g., 
IWim] [20091) . 

Stellar activity adds yet another dimension of complex- 
ity to the problem, because a (potentially evolving) sur- 
face brightness distribution also affects the transit profiles. 
The local brightness on the surface can either be decreased 
by dark spots or increased by bright faculae compared to 
the undisturbed photosphere. Spots (or faculae) located 
within the eclipsed section of the stellar surface lead to a 
decrease (increase) in the transit depth, and the actual pro- 
file depends on the distribution of those structures across 
the planetary path. Spots and faculae located on the non- 
eclipsed section of the surface do not directly affect the 
transit profile but change the overall level of the lightcurve. 
As transit lightcurves are, however, usually normalized with 
respect to the count rate level immediately before and af- 
ter the transit, the non-eclipsed spot contribution enters (or 
can enter) the resulting curve as a time-dependent modu- 
lation of the normalized transit depth. 

3.2. Transit lightcurve normalization 

As mentioned above, the normalization may affect the 
shape of the transit profiles. We now discuss two normal- 
ization approaches and compare their effect on the transit 
profiles. Let fi be the measured flux in time bin i, rii an 
estimate of the count rate level without the transit (hence- 
forth referred to as the 'local continuum'), and p a mea- 
sure of the unspotted photospheric level in the lightcurve, 
i.e., the count rate obtained in the respective band, when 
the star shows a purely photospheric surface. Usually, the 
quantity 

Vt = /i/n-i (2) 

is referred to as the 'normalized flux'. 

Normalization according to Eq.[2|can, however, result in 
variations of the transit lightcurve depth in response to non- 
uniform surface flux distributions as encountered on active 
stars. Assume a planet transits its host star twice. During 
the flrst transit the stellar surface remains free of spots, 
but during the second transit there is a large active region 
anywhere on the star not covered by the planetary disk (but 
visible). Consequently, the local continuum estimate, n^, 
for the second transit is lower, and the normalized transit 
appears deeper, although it is exactly the same transit in 
absolute (non-normalized) numbers. 

To overcome this shortcoming, we define the alternative 
normalization 

Zi = hi. (3) 

P 

In both cases the transit lightcurve is normalized with re- 
spect to the local continuum either by division or by sub- 
traction. The conceptual difference lies in the treatment of 
the local continuum level and how it enters the normal- 
ized transit lightcurve. Using Eq. [3]the observed transit is 
shifted, normalized by a constant, and shifted again. While 
the scaling in this case remains the same for all transits, the 
scaling applied in Eq.[5|is a function of the local continuum. 

Following the above example, we assume the same tran- 
sit is normalized using both Eq. [2| and Eq. [3] To evaluate 
the differences between the approaches we consider the ex- 
pression 

^ ift-ni)/p+l ^ ^ 
Vi fi/rii ~ 
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For Ui = p Eq. 3] holds as a strict equality, i.e. both nor- 
malizations yield identical results. The inequality equates 
to true if p > Ui and rii > fi - The first condition reflects the 
fact that the local continuum estimate should not exceed 
the photospheric lightcurve level, and the second one says 
that the lightcurve level is below the local continuum. The 
second condition is naturally fulfilled during a transit, and 
the first one is also met as long as faculae do not domi- 
nate over the dark spots during the tr ansit. At least in the 
case of CoRoT-2a lLanza et al.l ()2009() find no evidence for 
a significant flux contribution due to faculae, so that we 
conclude that the normalized transit obtained using Eq. [3] 
is always shallower than the one resulting from Eq. [2] unless 
Ui = p, in which case the outcomes are equal. 



3.2.1. Quantifying the normalization induced difference in 
transit depth 

Let us now study a single transit and consider data points 
covered by index set j, for which the term nj — fj reaches 
a maximal value of Tq at some index value j ~ T. At this 
position the normalization obtained from Eq. [3] is given by 
zt = (/t — nT)/p -I- 1 = —Tq/p + 1, whereas Eq. [2] yields 
Ut = JtI^t — {riT — Tq) / riT ■ These values are now used to 
compare the transit depths provided by the two normaliza- 
tions. Note that we assume here that the normalized depth 
is maximal at index T; this is always true for Eq. [3]but not 
necessarily for Eq. [5] which we regard a minor issue. Again 
we find zt = yr if — p. If, however, the local continuum 
estimate is given by ht ~ ap {a < 1), the results differ by 



ZT - IJT = Top 



(a- - 1) 



(5) 



Using the extreme values observed for CoRoT-2a {a ~ 
0.96, and Tq « 0.03 x p), the right hand side of Eq.lHyields 
« 1.3 X 10"'^ for the difference in transit depth, caused ex- 
clusively by applying two different normalization prescrip- 
tions. 



3.2.2. Which normalization should be used? 

For planetary research it is important to 'clean' the transit 
lightcurves of stellar activity in order to find the 'undis- 
turbed' profile associated with the planet only. As tran- 
sit lightcurves normalized using Eq. O are all scaled using 
the same factor, they preserve their shape and depth (at 
least relative to each other) and can, therefore, be com- 
bined consistently, which is not necessarily the case when 
Eq. [5] is used. Note that this does not mean that the ob- 
tained transit depth is necessarily the 'true' depth, because 
Eq.Oincludes the photospheric brightness level, p, as a time 
independent scaling factor. At least in the context of the 
lightcurve analysis, p cannot be determined with certainty 
as the star may not show an undisturbed surface during the 
observation; it may actually never show it. 

A problem evident in CoRoT lightcurve analyses is the 
existence of long term instrumental gradients in the data 
(cf. Sect. [2]). Modehng these trends by a 'shding' response, 
Rd, of the detector so that the relation between 'true' pho- 
tometry, Ci, and observation, Co.i, is given by Ci^o — Ci ■ Rd.i, 



Eq. [2 yields 



Obviously, the true photometry is recovered when the em- 
braced term equates to one. Yet the scaling of c in Eq. [1] is 
arbitrary, so that this is not necessarily the case. As long as 
qi, however, appropriately represents the shape of Rd.i the 
term provides a global scaling which cancels out in both 
Eq.Hand Eq. H 

For our transit analysis we argue in favor of the nor- 
malization along Eq. [31 We estimate the photospheric level 
from the highest count rate during the most prominent 
global maximum (at JD« 2454373.3) in each individual 
band. These estimates are based on the reduced lightcurves, 
in particular, the instrumental trend and the stellar com- 
panion have been accounted for. Throughout our analy- 
sis we use the values Pwhue = 703 000, pred = 489 000, 
Pgreen = 88 500, and Pblue = 124 500 (in units of e~/32s). 
Since even at that time spots are likely to be present on the 
stellar disk, these estimates might in fact represent lower 
limits to the true value of p. 

3.3. Transit profiles in CoRoT-2a 

The global lightcurve of CoRoT-2a shows pronounced max- 
ima and minima an d a temporally varia ble amplitude of the 
global modulation ([Alonso et al. [20081 ). It is natural to ex- 
pect the spot coverage on the eclipsed section of the stellar 
surface to be smallest where the global lightcurve is found 
at a high level, and transit events occurring during those 
phases should, thus, be least contaminated with the effects 
of stellar activity. The contrary should be the case for tran- 
sits during low lightcurve levels. 

To quantify the impact of activity on the transit profile, 
we define the transit equivalent width (TEW) of transit n 



TEWn 



(l-z„(t))dt«^(l-z„,,) 



St, 



(7) 



Ccorr.i — ■ I Rd.i 

It 



(6) 



where ti and tiv must be chosen so that they enclose the 
entire transit. Note that extending the integration bound- 
aries beyond the actual extent of the transit does not 
change the expectation value of Eq. [71 but only introduces 
an extra amount of error. The nominal unit of the TEW is 
time. 

3.3.1. A relation between transit equivalent width and global 
lightcurve modulation 

As outlined above we expect a larger impact of activity 
where the overall lightcurve level is low. When this is true, it 
should be reflected by a relation between the transit equiv- 
alent width and the transit continuum level (the overall 
lightcurve level at transit time). 

In Fig. m we show the distribution of TEWs as a func- 
tion of the local continuum level for all 79 transits ob- 
served with a 32 s sampling. There is a clear tendency 
for larger TEWs to be associated with higher continuum 
levels, thus, providing obvious evidence for activity-shaped 
transit lightcurves. In the same figure, we also show the 
best fit for a linear model relation, which has a gradient of 
d(TEW)/d(CL)= (5± 1.5) x 10"'' s/(e-32 s). 

To corroborate the reality of the above stated correla- 
tion, we calculated the correlation coefficient, R. Its value of 
R = 0.642 confirms the visual impression of a large scatter 
in the distribution of data points (cf., Fig.[T]). We estimate 
the statistical error for a single data point to be « 0.1%, so 
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Mean continuum level [e7(32 s)] 



Fig. 1. Transit equivalent width (TEW) vs. transit contin- 
uum level as well as the best-fit linear model. 



that the scatter cannot be explained by measurement er- 
rors. To check whether the continuum level and the TEW 
are independent variables, we employ a t-test and find the 
null hypothesis (independent quantities) to be rejected with 
an error probability of 1.8 x 10"^", so that the correlation 
between the TEWs and the continuum level must be re- 
garded as highly significant. 

As a crosscheck for the interpretation of this finding, 
we also investigated the distribution of TEWs against time, 
which shows no such linear relation {R ~ 0.110). Therefore, 
we argue that the effect is not instrumental or caused by 
our data reduction, but physical. 



3.4. Comparing high and low continuum level transits 

Since activity is evident in the profiles of the transit 
lightcurves, we further investigate its effect by comparing 
the most and least affected transit lightcurves. Therefore, 
we average the ten transits with the highest continuum lev- 
els (no. 3, 16, 42, 47, 50, 55, 68, 73, 76, and 81) and compare 
the result to an average of the ten transits with the low- 
est continuum level (no. 15, 23, 35, 40, 43, 69, 72, 75, 77, 
and 80). In Fig. [2] we show the two averages as well as our 
computed lower envelope (see Sect. 13.51) superimposed on 
the entire set of folded photometry data points. The dis- 
tribution of the entire set is denoted by a color gradient 
(red) with stronger color indicating a stronger concentra- 
tion of data points. The curve obtained from the transits 
at 'low continuum state' is clearly shallower, as was already 
indicated by the TEW distribution presented in Fig. [TJ 

The difference in TEW amounts to « 15.5 s in this 
extreme case. We checked the significance of this number 
with a Monte Carlo approach. On the basis of 20 randomly 
chosen transits, we constructed two averaged lightcurves 
using 10 transits for each and calculated the difference in 
TEW. Among 1 000 trials we did not find a single pair with 
a difference beyond 12 s, so that the result is likely not due 
to an accidental coincidence. 



3.5. Obtaining a lower envelope to the transit profiles 

As was demonstrated in the preceding section, activity 
shapes the transit lightcurves, and we cannot exclude that 
every transit is affected so that a priori no individual pro- 
file can be used as a template representing the 'undisturbed' 
lightcurve. The distortion of the individual profiles is, how- 
ever, not completely random, but the sign of the induced 
deviation is known as long as we assume that dark struc- 
tures dominate over bright faculae which seems justified 



for CoRoT-2a (jLanza et al.ll2009f) . In this case activity al- 
ways tends to raise the lightcurve level and, thus, decreases 
the transit depth. Therefore, the best model for the undis- 
turbed profile can be obtained as a lower envelope to the 
observed transit profiles. 

Assume a set of Nt transit observations, with the as- 
sociated photometry folded at a single transit interval pro- 
viding the set LCx^i of transit data points. If the lower en- 
velope were observed (which might even be true) , it would 
in principle look like every other lightcurve. In particular, 
it shows the same amount of intrinsic scattering (not in- 
cluding activity), characterized by the variance (Tq. 

Suppose we have an estimate of the variance 



1 ^ 

^0 ~ j;^Y.^lCt,j - 



(8) 



with the (unknown) expectation value iij and the number 
of data points N. The aim of the following effort is to iden- 
tify the lowest conceivable curve sharing the same variance. 
In order to achieve this, we divide the transit span into a 
number of subintervals each containing a subsample, s, of 
LCt- The distribution of data points in s is now approxi- 
mated by a 'local model', lm{'y), with a free normalization 
7; Im can for instance be a constant or a gradient. Given Im, 
we adapt the normalization in order to solve the equation 



j:^H{lm{^)-LCs) 



(9) 



where H denotes the Heaviside function {H{x) — 1 for 
a; > and H{x) = otherwise). In this way we search 
for the local model compatible with the known variance 
of the lower envelope. The ratio on the left hand side of 
Eq. [5] represents a variance estimator exclusively based on 
data points below the local model. It increases (strictly) 
monotonically except for the values of 7 where the local 
model 'crosses' a data point and the denominator increases 
by one instantaneously. Therefore, there may be more than 
one solution to Eq. [9l From the mathematical point of view 
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Fig. 2. Average transit lightcurves obtained by combining 
the ten profiles which exhibit the highest (thick dashes) and 
lowest (thin dashes) continuum levels. The crosses indicate 
our lower envelope estimate and the color gradient (red) 
illustrates the distribution of data points for all available 
transits. 
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all solutions are equivalent, but for a conservative estimate 
of the lower envelope the largest one should be used. 

In Fig. [2] we show the lower envelope which is in much 
better agreement with the average of the high continuum 
transit profiles than with its low continuum counterpart. 
The derivation of the lower envelope is based on Eq. [9l 
To obtain an estimate of Cg we fitted a 500 s long span 
within the transit flanks (3 500 ± 250 s from the transit 
center) where activity has little effect with a straight line 
and calculated the variance with respect to this model. The 
resulting value (using normalized flux) of CTq = 1.6 x 10"^ 
was adopted in the calculation. Furthermore, we chose a 
bin width of 150 s, and the 'local model' was defined as a 
regression line within a ±100 s span around the bin center. 
Additionally, we postulated that at least 8 (out of « 350) 
data points per bin should be located below the envelope, 
which proved to make the method more stable against the 
effect of outliers and has little impact otherwise. 

3.6. Transit profiles in different color channels 

CoRoT observes in three different bands termed 'red', 
'green', and 'blue'. In the following we present a qualita- 
tive analysis of the transit profiles in the separate bands. 
In the case of CoRoT-2a approximately 70% of the fiux is 
observed in the red band, and the remaining 30% is more 
or less equally distributed among the green and blue chan- 
nels. In order to compare the profiles we average all avail- 
able transits in each band individually and normalize the 
results with respect to their TEW, i.e., after this step they 
all have the same TEW. The resulting profiles represent the 
curves which would be obtained if the stellar fiux integrated 
along the planetary path was the same in all bands. 

In Fig. [3] we show the thus normalized transit lightcurves 
(TEW=1) obtained in the three bands. The normalized 
transits show a difference both in their fiank profile and 
their depth. The blue and green transit profiles are nar- 
rower than the red one, and they are deeper in the center. 
This behavior is most pronounced in the blue band so that 
the green transit lightcurve virtually always lies in between 
the curves obtained in red and blue. 
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Fig. 3. Left panel: Normalized transit in the three CoRoT 
bands red, green, and blue obtained by averaging all avail- 
able data. Upper right: Close-up of the transit center. Lower 
right: Close-up of the ingress fiank of the transit. 



The behavior described above can be explained by a 
color-dependent limb darkening law, with stronger limb 
darkening at shorter wavelengths as predicted by atmo- 
spheric models (1 Claret 2 0041 ) and observed on the Sun 
(jPierce fc Slaughter ,1977j). Note t hat we checked that an- 
alytical transit models ( |Pa]||2008D generated for a set of 
limb darkening coefficients, indeed, reproduce the observed 
behavior when normalized with respect to their TEW. 

Normalizing the averaged transits not with respect to 
TEW but using Eq. [3] yields approximately the same depth 
in all bands, while the difference in the fianks becomes more 
pronounced. The reason for this could be a wrong relative 
normalization, which can e.g. occur if the eclipsed section 
of the star is (on average) redder than the rest of the sur- 
face due to pronounced activity or gravity darkening, or 
it may be a relic of an inappropriate treatment the com- 
panion's flux contribution. Whatever the explanation, it is 
clear from Fig. [3] that the flanks and the centers in the 
individual bands cannot be reconciled simultaneously by 
a re-normalization. Therefore, our analysis shows that the 
transit lightcurves are color dependent. 

4. Stellar activity and planetary parameters 

The preceding discussion shows that stellar activity does 
have a considerable influence on the profile of the transit 
lightcurves, and, therefore, the derivation of the planetary 
parameters will also be affected. Below we determine the 
radius and the orbit inclination of CoRoT-2b taking activ- 
ity into account and discuss remaining uncertainties in the 
modeling. 

4.1. Deriving the planetary radius and inclination from the 
lower envelope profile 

In the analysis presented bv I Alons"o" et al.l daOOS"), the fit of 
the planetary parameters is based on the average of 78 tran- 
sit lightcurves (see Table [1] for an excerpt of their results). 
While this yields a good approximation, the results still in- 
clude a contribution of stellar activity, and an undisturbed 
transit is needed to calculate 'clean' planetary parameters. 

We follow a simplified approach to estimate the impact 
of activity on the planetary parameters. In particular, we 
use the lower envelope derived in Sect. 13.41 as the best avail- 
able model for the undisturbed t r ansit . Starting from the 
results reported bv lAlonso et al.l (|2008l ). we re-iterate the 
fit of the planetary parameters. In our approach we fix the 
parameters of transit timing, the semi-major axis and stel- 
lar radius , and the limb darkening coefficients at the values 
given bv lAlonso eFall (|2008, ) (cf. Table [J). The two free 
parameters are the planetary radius and its inclination. 

Note that limb darkening coefficients recovered by 
lightcurve analyses are not relia ble especiall y when more 
than one coefRcient is fitted (e.g., Winn 2009). However, as 
an accurate calibration of the CoRoT color ba nds is not yet 
availa ble and the coefficients determined by lAlonso et al.l 
(120081) roughly correspond to numbers predic ted by stellar 
atmosphere models Q, we decided to use the lAlonso et al.l 
values here, too. This, furthermore, simplifies the compar- 
ison of the results. 

^ For T eff = 5 6 00 K and log(g) = 4.5 the PHOENIX models 
given by IClareti (|2004f ) yield quadratic limb darkening coeffi- 
cients of Ua ~ 0.46 and ut = 0.25 in the Sloan-r' band. 



6 



S. Czesla et al.: How stellar activity affects the size estimates of extrasolar planets 



Fo r the fit, we use the analytical models given by iPall 
(|2008f ) in combination with a Nelder-Mead simplex algo- 
rithm fe.g.. [Press,199a . 
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Fig. 4. Lower envelope of all normalized transit lightcurves 
(already shown in Fig. [2) and our model fit. 

The result of our modeling is illustrated in Fig. ID The 
best radius ratio is Rp/Rg = 0.172 ±0.001 at an inclination 
of 87.7°±0.2°. Note that the quoted errors are statistical er- 
rors, only valid in the context of the model. These numbers 
need to be compared to the values Rp/Rs = 0.1667±0.0006 
and 87.84°±0.1° (cf., Table[T]) derived without taking activ- 
ity effects into account. The best fit inclination is com pati- 
ble with the value determined bv lAlonso et al.l ()2008f ). but 
'our' planet is larger by « 3%. The planet's size mainly 
depends on the transit depth which is, indeed, affected at 
about this level by both normalization (Sect. 13.2. ip and 
stellar activity (Sect. 1?^ . 

Clearly, the derived change in Rp/Rg of 0.005 is much 
larger than the statistical error obtained from light curve 
fitting, and, therefore, the neglect of activity leads to sys- 
tematic errors in excess of statistical errors. While the over- 
all effect in planet radius is « 3%, the error in density be- 
comes ~ 10%. Such errors are certainly tolerable for mod- 
eling planetary mass-radius relationships, they are not tol- 
erable for precision measurements of possible orbit changes 
of such systems. 

4.1.1. Planetary parameters and photospheric level 

As already indicated the normalization according to Eq. [3] 
relies on a 'photospheric lightcurve level', p, which enters 
it as a global scaling factor and, therefore, also interferes 
with pinning down the planet's properties. 

In a simple case the star appears as a sphere with 
a purely photospheric surface, and the observed transit 
depth, /o, can be identified with the squared ratio of plan- 
etary and stellar radii 

^ _ max(n, ^ /,) ( RpV r /^n^ 

Here Ld is a correction factor accounting for limb darken- 
ing. Yet, when the observed star is active and the lightcurve 
is variable, there is no guarantee that the maximum point in 
the observed photometry is an appropriate representation 
of the photospheric stellar luminosity. Persistent inhomo- 
geneities such as polar spots and long-lived spot contribu- 
tions modulate the lightcurve, so that the pure photosphere 
might only be visible anytime the star is not observed or 
possibly never. 



Assume our estimate, Pmi of the photospheric level un- 
derestimates the 'real' value, p, by a factor of < c < 1 so 
that Pm = p ■ c and fie^i denotes the lower envelope transit 
lightcurve. Then the measured transit depth, /„, becomes 

J. _ max(nt - fie^) _ / Rp\^ Ld ^^^^ 
™ Pm \R* J c 

and another scaling factor must be applied to the radius 
ratio. While p„i is a measured quantity, c is unknown, and 
if we neglect it in the physical interpretation, i.e., the right 
hand side of Eq. [11] the ratio of planetary and stellar radii 
will be overestimated by a factor of 1/ y/c. 

The value of c cannot be quantified in the frame of this 
work; only an estimate can be provided. Doppler imaging 
studies revealed that polar spots are common, persisten t 
structures in young, active stars (e.g. iHuber et al.ll2009[ ). 
Assuming that polar spots also exist on Corot-2a and that 
they reach down to a latitude of 70°, they occupy roughly 
2% of the visible stellar disk. Adopting a spot contrast of 
50%, c becomes 0.99 in this case, and the planet size would 
be overestimated by 0.5%. As the poles of CoRoT-2a are 
seen under a large viewing angle, their impact would, thus, 
be appreciably smaller than the amplitude of the global 
brightness modulation (ca. 4%). Nonetheless, in sign, this 
effect counteracts the transit depth decrease caused by ac- 
tivity, and if the polar spots are larger or symmetric struc- 
tures at lower latitudes contribute, it may even balance it. 

5. Discussion and conclusion 

Stellar activity is clearly seen in the CoRoT measured tran- 
sit lightcurves of CoRoT-2a, and an appropriate normal- 
ization is necessary to preserve as much of the true transit 
lightcurve profile as possible. 

The transit profiles observed in CoRoT-2a are affected 
by activity, as is obvious in many transits where ac- 
tive regions cause distinct 'bumps' in the lightcurve (e.g. 
IWolter et al]l2009l) . Yet, our analysis reveals that not only 
profiles with bumps but presumably all transit profiles are 
infiuenced by stellar activity. This is evident in a relation- 
ship between the transit equivalent width and the level 
of the global lightcurve: Transits observed during periods 
where the star appears relatively bright are deeper than 
those observed during faint phases. We demonstrated this 
correlation to be extremely significant, but also the data 
points show a large scatter around an assumed linear model 
relation. If the star were to modulate its surface bright- 
ness globally and homogeneously, this relation would be 
perfectly linear except for measurement errors. Therefore, 
we interpret the observed scatter as a consequence of sur- 
face evolution. When the global lightcurve is minimal, we 
also find more spots under the eclipsed portion of the sur- 
face, but only on average and, for an individual transit, this 
needs not to be the case. Thus, the surface configuration is 
clearly not the same for every minimum observed. 

Additionally, we demonstrated that the transit profiles 
show a color dependence compatible with a color-dependent 
limb darkening law as expected from stellar atmospheric 
models and the solar analogy. 

All these influences potentially interfere with the de- 
termination of the planetary parameters. Using our lower 
(white light) transit envelope, we determined new values 
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for the planet-to-star radius ratio and the orbital inclina- 
tion. While the latter remains compatible with previously 
reported results, the planet radius turns out to be larger 
(compared to the star) by about 3%. Although our ap- 
proach accounts for many activity related effects, a num- 
ber of uncertainties remain. For example the photospheric 
lightcurve level needed for transit normalization cannot be 
determined with certainty from our analysis and the same 
applies to the limb darkening law. More certain than the 
planetary parameters themselves is, therefore, our conclu- 
sion that the errors in their determination are much larger 
than the statistical ones. 

While CoRoT-2a is certainly an extreme example of an 
active star, stellar activity is a common phenomenon espe- 
cially on young stars. Therefore, in general, stellar activity 
cannot be neglected in planetary research if the accuracy 
of the results is to exceed the percent level. 
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